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We study the renormalization of the effective mass and trigonal warping of bilayer graphene by the 
electron- electron interaction. One consequence of such a renormalization in the low-energy bands 
of a bilayer crystal consists of a small reduction of the critical density of the Lifshitz transition (the 
crossover between the single-pocket and four-pocket topology of the Fermi surface). 

PACS numbers: 71.10.Ca, 71.18.+y, 71.70.Gm 



The bilayer graphene crystal is one of several allotropic 
forms of carbon one can fabricate using micro-mechanical 
cleavage of bulk graphite^. The crystalline structure of 
bilayer graphene is derived from the Bernal stacking 2 of 
layers in the "mother" graphite crystal. It is shown in 
Fig. H{a) as two hexagonal carbon lattices with inequiv- 
alent sites A, B (bottom layer) and A, B (top layer) such 
that the A sites are above the B sites, while A (B) is 
located above (below) the midpoint of a hexagon on the 
other layer. Theoretical studies 3 , transport 1 , ARPES 4 
and optical 5 characterization of bilayer graphene have 
indicated that this material is a gapless semiconductor 
with two split-bands separated by ±71, where 71 is the 
closest-neighbor interlayer (AB) coupling, and two "low- 
energy" almost parabolic bands with effective mass de- 
termined by both 71 and the intralayer coupling 70. 

Detailed tight-binding model studies 3 of bilayer 
graphene taking into account the next-neighbor (B A) in- 
terlayer coupling 73 have indicated that the dispersion 
of the low-energy bands in it is strongly anisotropic at 
small momenta p ~ 7371 /(7o&), with the anisotropy pa- 
rameterized by ^3 = -^ItP". As a result, the Fermi line 
in bilayer graphene may undergo a topological (Lifshitz) 
transition 6 : from singly connected at high carrier den- 
sity, n e fh) > n L to four separate pockets in momentum 
space, for n e (h) < n L- It is common for two-dimensional 
electron systems^ that electron-electron (e-e) repulsion 
renormalizes the single-particle dispersion of carriers. In 
monolayer graphene, where electrons have a characteris- 
tic Dirac spectrum, this leads 8 to an increase of the Dirac 

velocity, from v = ^fp to v(p) = v (l + f ln(A/p)) 
with a = e 2 /e s v. In this Communication we investigate 
the effect of the e-e interaction on the split-band gap 71 , 
the effective mass m, and the dispersion anisotropy pa- 
rameter vs. We show for the bilayer that in the Hartree- 
Fock theory the e-e repulsion increases 71 and V3, but 
renormalizes the effective mass 9 and the Lifshitz transi- 
tion density downwards. Also, we find that, in contrast 
to the monolayer, these corrections to the electronic spec- 
trum do not contain infrared divergencies, and that the 
renormalization of ul is weak, so that the single-particle 
tight-binding model gives a realistic estimation^ for the 
Lifshitz transition density in bilayers. 
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FIG. 1: (Color online) (a) The crystalline structure of bilayer 
graphene. (b) The first Brillouin zone and the constant energy 
cuts about the K, K' points. The dark shaded area represents 
the four pockets in momentum space filled by electrons in the 
gas with \n\ < ul below the Lifshitz transition density, the 
light shaded area is the simply connected Fermi line for higher 
densities. 



The theory presented below is based on the four-band 
model 3 - describing electrons in bilayer graphene in the 
vicinity of the valley centers corresponding to the K 
(£ = 1) and K' (£ = —1) first Brillouin zone corners 
[see Fig. [U^b)]. The 4x4 tight-binding Hamiltonian, 
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which will be the subject of renormalization by the e- 
e interaction, is written in the basis of sublattice Bloch 
states [ipAi ^S'^Ai ^] in valle y K and fe 5 ^a, fe, 
in valley K r . Notice that H has a natural ultraviolet 
momentum cutoff about A ~ 1/a. This Hamiltonian de- 
termines four bands 3 , ±ei 5 2(p). Two split-off bands ±62 
start at energies ±71. For small values of momentum, 
vp <C 7i/4, the two low-energy bands dzei that touch 
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each other at Kj K' can be attributed to a 2 x 2 effective 
Hamiltonian 3 

where m = ^2, and H' acts on [i/jai ^g] i n valley if and 
^a] in valley if 7 . The second term in H r causes a tri- 
angular distortion of the electronic dispersion illustrated 
in Fig. QJb), where the Lifshitz transition of the electron 
Fermi line topology is explained, as a singly-connected 

line ei(p) = E splits into four disconnected pockets at 

2 

E < cl = i^r-, which would occur at a critical density 10 
(h restored) 

In the Hartree-Fock approximation, the change in the 
single-particle Hamiltonian H = H + E for electrons can 
be described using the self-energy diagram 

O 

£ = j + /'vfqfN (4) 



where the solid line denotes the causal Green's function of 
electrons in the bilayer and the dashed line stands for the 
interaction V(q) = 27re 2 /e s qx(q), where e s is the dielec- 
tric constant of the environment (e s = 1+ 2 4 ' 5 = 2.75 on a 
Si02 substrate). The first (Hartree) diagram in Eq. (|4]), 
together with the positive background charge cancels the 
momentum-independent energy shift determined by the 
second (Fock) term. Below, the contribution from the 
Fock term is estimated both for the bare Coulomb inter- 
action, V(q) = 27re 2 /e s q and for the screened Coulomb 
interaction, where the dielectric function x(q) ls evalu- 
ated in the random phase approximation (RPA). The 
self-energy E is evaluated for bilayer graphene at zero 
doping. Since the Lifshitz transition occurs at very low 
doping level, using the renormalized value of the bilayer 
parameters that result from this approximation is justi- 
fied as long as the critical Fermi energy of this transition 
remains small in comparison to the energy scale 71. 



In the 4x4 representation the self-energy is a matrix, 
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Here, we omitted two unimportant terms, and define n = q x + iq y , ir = p x + ip y = pe*^, and ei 5 2 = ^1,2 (q + p) with 
4(6 P) = y + (« 2 + y) P 2 + (- l T\J\ (7? " « 3 V) 2 + «V (7? + « 3 V) + 2£ 7lW3 «V cos3^. (6) 



The self-energy E in Eq. ([5]) is ^-independent, as usual 
in the Hartree-Fock approximation. Schematically, the 
result of integration over q in Eq. ([5]) can be represented 
as 



£(£,p)=e 2 
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where A = £A l7 r + 0(p 2 ), C = fC + 0(p 2 ), B = 
^B 1 tt + ^B 2 n 2 + 0(p 3 ), and Z = ^M l7 r + 0(p 3 ). 

Here we expanded E in powers of p, in order to keep only 
the terms relevant for the renormalization of the bilayer 



parameters v, ^3,7i,m and ul- Term A in Eq. takes 
into account the correction to the intra-layer velocity v. 
A similar effect is known to occur in monolayers 8 . Term 
C in Eq. (0) modifies the inter-layer coupling 71. The 
part of terms B and Z that is linear in 7r, 7f yields renor- 
malization of vs. The quadratic part of term B combines 
with A and C to determine the renormalization of the 
effective mass, 
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which is a result similar to that obtained by Borghi et al 9 . 
Together, all of these determine the shift of the Lifshitz 
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bare Coulomb (any e s ) 


screened suspended (e s = 1) 


screened on Si02 (e s = 2.75) 




(0.58-^0.66) + f In (vA/71) 


(0.03 0.04) + 0.032 In (vA/71) 


(0.076 -r- 0.098) + 0.071 In (vA/71) 


B 1 


-0.063 


-0.0013 


-0.0034 


B 2 


0.38 


0.016 


0.039 


Co 


(0^0.36) + | In (vA/71) 


(-0.025 -r- -0.14) + 0.032 In (vA/71) 


(-0.047 -r- -0.021) + 0.071 In (vA/71) 


Mi 


(0.7-=-l.l) + iln(i;A/7i) 


(0.036 -r- 0.082) + 0.032 In (vA/71) 


(0.09 0.19) + 0.071 In (vA/71) 



TABLE I: The physical parameters of the self-energy £(£,p) (Eq. J7}), for the bare and screened Coulomb interaction. Where 
two limits are given, they correspond to integrating the finite (not logarithmically divergent) part up to q = ^ and to q = 00, 
respectively. 



transition density, 

n L = n L (1 + 2a (C + Mi + B ± - 2A ± )) . (8) 

The numerical values of the coefficients Ai, B2, 
Co, and M\ in Eq. ([7j) were calculated for suspended 
graphene flakes and flakes on Si02 substrate with and 
without screening of the Coulomb interaction taken into 
account^ (Table H}. 

Screening of the e-e interaction in graphene is taken 
into account in the RPA, by the effective static dielectric 
function 

x( q ) =1 - nq)n(q), 

H(q)=-< / ^Tv(G (p+|, W )G (p-f, W )), 

where II (q) is the polarizability of the electron gas. Fig. [2] 
shows the numerical evaluation of 11(g), which deviates 
from the constant value 12 ~ — 0.447i/v 2 obtained earlier 
using the two-band Hamiltonian H' at q ~ Ji/2v, where 
the band has a crossover from parabolic to an almost 
linear behavior. To obtain analytical asymptotic expres- 
sions in the relevant limits first, consider q > q* > ji/v, 
where the dominant contribution comes from p ^> 71 /v, 
and the free electron Green's function reduces to the free 
Green's function of two decoupled monolayers. Then 8 , 
n = -±Nq/v with N = 8, and 
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For q < q*, where the electron dispersion is mostly 
parabolic, screening has the form of Thomas- Fermi 
screening with the radius tq = 0.1 7v/ji (for e s = 1). 
A sufficient simultaneous analytical description of both 
regimes can be done using an interpolation formula, 



V(q<q*) = ((3 + 0q 2 +r 1 q A )e- 



(10) 



where the parameters ^,7,7^,^ are determined by a nu- 
merical fit shown in Fig. [2j and the value of q* is deter- 
mined from the intersection of V(q < q*) and V(q > q*). 
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FIG. 2: (Color online) The renormalized interaction 11 
V(q) from numerics (symbols), its large momentum asymp- 
totics Eq. (|9} (dashed line) and low-momentum asymptotics 
Eq. (|1Q|) (solid line) for e s = 1 and 2.75. Inset: the static 
polarization H(q) from numerics and its asymptotics. Table: 
parameters of the small q approximate interaction (Eq. (|10[0. 



Finally, we determine that the bilayer parameters are 
renormalized by the e-e repulsion as 
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where the numerical values of all J's are listed in Ta- 
ble [III This result shows that the intralayer velocity 

the interlayer hopping 71, and the trigonal distortion 
parameter ^3 are all enhanced by the e-e interaction^. 
With bare Coulomb interaction this enhancement would 
be grossly overestimated if we used the actual value 
a = e 2 /hv w 2.19 for graphene. 

However, after screening is taken into account, we find 
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a 


S v 


S V3 


K 




-l 


I 


hve s 


0.58^ 


0.66 


0.67^ 


1.05 


0^0.36 


1.00-=- 


1.51 


II 


0.28 


0.07^ 


0.09 


0.08^ 


0.18 


±0.05 


0.17 -r- 


0.20 


III 


0.20 


0.06^ 


0.08 


0.07^ 


0.15 


±0.04 


0.17 -r- 


0.19 



TABLE II: Constants for Eq. ([TT]) for bare Coulomb interac- 
tion (I), screened Coulomb interaction in suspended graphene 
(II), and screened Coulomb interaction in graphene on SiC^ 
substrate (III). 11 Where two limits are given, they correspond 
to integrating the finite (not logarithmically divergent) part 
up to q — ji/v and to q = oo, respectively. 
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FIG. 3: (Color online) The relative reduction Sul/til — 
ul/til — 1 of the Lifshitz transition density by the e-e in- 
teraction for a suspended sheet and one on SiCb substrate. 
The vertical line shows the estimation for A « 70 /v. 



that one has to usein Eq. (fTTj) 5 ~ 0.28 for a suspended 
flake and a ~ 0.20 for a flake on SiC>2, which gives a much 
weaker renormalization effect than an estimate using bare 



Coulomb interaction. Also, using Eqs. ([8]) and ([TT]h one 
may see that in the theory ignoring screening the renor- 
malization of the Lifshitz transition density ul would 
be hugely overestimated, resulting in its disappearance 
(ul — > 0). However, having taken into account the effect 
of the reduction of the e-e repulsion by screening, we find 
a much smaller shift in the value of til from that deter- 
mined using the tight-binding model for noninteracting 
electrons^^. The calculated shift in the transition den- 
sity til for the screened Coulomb interaction is shown 
in Fig. [3] as a function of momentum cutoff A, with the 
vertical line corresponding to A « 70/^. This determines 
an approximately 15% reduction of the Lifshitz transi- 
tion density for a suspended sheet and about 12% for a 
bilayer on a Si02 substrate. 

The result of the above-presented analysis of the bi- 
layer band parameters and the Lifshitz transition density 
til suggests that their renormalization by the e-e repul- 
sion is relatively weak, due to the screening of the e-e 
interaction by the electrons themselves. Thus, we con- 
clude that the Lifshitz transition is not trivially hindered 
by many-body effects. The Lifshitz transition can be de- 
tected through the singularity of the thermopower— that 
develops when a neck forms from the central pocket to 
the side pockets of the Fermi line in bilayer graphene 
(Fig. DIb)). Such an observation would require samples 
of high homogeneity, but in contrast to bulk metals for 
which the Lifshitz transition was first discussed, such an 
experiment would be possible in graphene since the car- 
rier density in graphene can be directly controlled using 
external gates. 
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